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Short Version 
Abstract 

We show that a PL-reahzation of a closed connected manifold of dimension n — 1 
in M" (n > 3) is the boundary of a convex polyhedron if and only if the interior of each 
(n — 3)-face has a point, which has a neighborhood lying on the boundary of a convex 
n-dimensional body. This result is derived from a generalization of Van Heijenoort's 
theorem on locally convex manifolds to the spherical case. No initial assumptions about 
the topology or orientability of the input surface are made. Our convexity criterion for 
PL-manifolds implies an easy polynomial-time algorithm for checking convexity of a 
given PL-surface in M". The algorithm is worst case optimal with respect to both the 
number of operations and the algebraic degree. The algorithm works under significantly 
weaker assumptions and is easier to implement than convexity verification algorithms 
suggested by Mehlhorn et al (1996-1999), and Devillers et al.(1998). A paradigm of 
approximate convexity is suggested and a simplified algorithm of smaller degree and 
complexity is suggested for approximate floating point convexity verification. 

There is a number of theorems that infer global convexity from local convexity. The 
oldest one belongs to Jacque Hadamard (1897) and asserts that any compact smooth surface 
embedded in M^, with strictly positive Gaussian curvature, is the boundary of a convex 
body. Local convexity can be defined in many different ways (see van Heijenoort (1952) 
for a survey). We will use Bouligand's (1932) notion of local convexity. In this definition a 
surface M in the affine space is called locally convex at point p if p has a neighborhood 
which lies on the boundary of a convex n-dimensional body Kp] if /^p\p lies in an open 
half-space defined by a hyperplane containing p, M is called strictly convex at p. 

This paper is mainly devoted to local convexity of piecewise-linear (PL) surfaces, in 
particular, polytopes. A PL-surface in R" is a pair M = (M, r), where M is a topological 
manifold with a fixed cell-partition and r is a continuous realization map from M to R" that 
satisfies the following conditions: 
1) r is a bijection on the closures of all cells of M 

1) for each fc-cell C of M the image r{C) lies on a /c-dimensional affine subspace of M"; r(C) 
is then called a fc-face of M. 
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Thus, r need not be an immersion, but its restriction to the closure of any cell of M 
must be. By a fixed cell-partition of M we mean that M has a structure of a CW-complcx 
where all gluing mappings are homeomorphisms (such complexes are called regular by J.H.C. 
Whitehead). All cells and faces are assumed to be open. We will also call M = (M, r) a 
PL-reahzation of M in R". 

Definition 1 We say that M = (M, r) is the boundary of a convex body P if r is a home- 
omorphism between M and dP. 

Hence, we exclude the cases when r(M) coincides with the boundary of a convex set, but 
r is not injective. Of course, the algorithmic and topological sides of this case are rather 
important for computational geometry and we will consider them in further works. Notice 
that for n > 2 a closed (n — l)-manifold M cannot be immersed into by a non-injective 
map r so that r(M) is the boundary of a convex set, since any covering space of a simply 
connected manifold must be simply connected. However, such immersions are possible in 
the hyperbohc space H". 

Our main theorem asserts that any closed PL-surface M immersed in R" (n > 3) with 
at least one point of strict convexity, and such that each {n — 3)-cell has a point at which 
M is locally convex, is convex. Notice that if the last condition holds for some point on an 
{n — 3)-face, it holds for all points of this face. 

This theorem implies a test for global convexity of PL-surfaces: check local convexity on 
each of the (n — 3)-faces. Notice that if for all k and every k-iace there is an (n— fc — l)-sphere, 
lying in a complementary subspacc and centered at some point of F, such that § fl F is a 
convex surface, then r is an immersion. The algorithm implicitly checks if a given realization 
is an immersion, and reports "not convex" if it is not. 

The pseudo-code for the algorithm is given in this article. The complexity of this test 
depends on the way the surface is given as input data. Assuming we are given the coordinates 
of the vertices and the poset of faces of dimensions n— 1, n — 2, n — 3 and 0, OR, the equations 
of the facets, and the poset of faces dimensions n — 1, n — 2, and n — 3, the complexity of 
the algorithm for a general closed PL-manifold is 0{fn-3,n-2) = 0{fn~3,n-i), where fk,i is 
the number of incidences between cells of dimension k and /. If the vertices of the manifold 
are assumed to be in a sufficiently general position, then the dimension of the space does 
not affect the complexity at all. Another advantage of this algorithm is that it consist of 
fn-3 independent subroutines corresponding to the {n — 3)-faces, each with complexity not 
exceeding O in the number of (n — l)-cells incident to the (n — 3)-facc. 

The complexity of our algorithm is asymptotically equal to the complexity of algorithms 
suggested by Devillers et al (1998) and Mehlhorn et al (1999) for simplicial 2-dimensional 
surfaces; for n > 3 our algorithm is asymptotically faster than theirs. These authors verify 
convexity not by checking it locally at {n — 3)-faces, but by different, rather global methods 
(their notion of local convexity is, in fact, a global notion). Devillers et al (1998) and 
Mehlhorn et al (1999) make much stronger initial assumptions about the input, such as the 
orientability of the input surface; they also presume that for each {n — l)-face of the surface 
an external nornial is given, and that the directions of these normals define an orientation 
of the surface. Then they call the surface locally convex at an (n — 2)-face F if the angle 
between the normals of two (n — l)-faces adjacent to F is obtuse. Of course this notion of 
"local" convexity is not local. 
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The main theorem is deduced from a direct generahzation of van Heijenoort's theorem to 
the spherical case. Van Heijenoort's theorem asserts that an immersion in M" of any closed 
connected manifold M, which is locally convex at all points, strictly locally convex in at 
least one point, and is complete with respect to the metric induced by r, is the boundary 
of a convex d-dimensional set. Van Heijenoort (1952) noticed that for n = 3 his theorem 
immidiately follows from four theorems contained in Alexandrov (1948); however, acoording 
to van Heijenoort, Alexandrov's methods do not extend to n > 3 and his approach is 
technically more complicated. We show that this theorem also holds for spheres, but not for 
the hyperbolic space. While all notions of affine convexity can be obviously generalized to 
the hyperbolic space, there are two possible generalizations in the spherical case; neither of 
these generalizations is perfect. The main question is whether we want all geodesies joining 
points of a convex set S to be contained in S, or at least one. In the first case subspheres are 
not convex, in the second case two convex sets can intersect by a non-convex set. The latter 
problem can be solved by requiring a convex set to be open, but this is not very convenient, 
since again, it excludes subspheres. We call a set in X" convex if for any two points p, q e S*, 
there is some geodesic [p, q] C 5. 

Proposition 2 // the intersection I of two convex sets in S", n > is not convex, then I 
contains two opposite points. 

To have unified terminology we will call subspheres subspaces. 

Besides the algorithmic implications, our generalization implies that any (n — 3)-simple 
PL-surface in with convex facets is the boundary of a convex polyhedron. 

1 van Heijenoort- Alexandrov's Theorem 
for Spaces of Constant Curvature 

Throughout the paper X" denotes M", S", or H". Following the original proof of van 
Hejeenoort's, we will now show that his theorem holds in a somewhat stronger form for 
S" for n > 2. van Hejeenoort's theorem does not hold for unbounded surfaces in H". We will 
give three different kinds of counterexamples and pose a conjecture about simply connected 
locally compact embeddings of manifolds in H". 

Imagine a "convex strip" in 3D which is bent in the form of handwritten cp so that it 
self-intersects itself, but not locally. Consider the intersection of this strip with a ball of 
appropriate radius so that the self-intersection of the strip happens to be inside the ball, 
and the boundary of the strip outside. Regarding the interior of the ball as Klein's model 
of we conclude that the constructed surface is strictly locally convex at all points and 
has a complete metric induced by the immersion into the hyperbolic space. This gives an 
example of an immersion of a simply connected manifold into which does not bound a 
convex surface. Notice that in this counterexample the surface self-intersects itself. 

Consider the (affine) product of a non-convex quadrilateral, lying inside a unit sphere 
centered at the origin, and a line in W\ The result is a non-convex polyhedral cylindrical 
surface P. Pick a point p inside the sphere, but outside the cylinder, whose vertical projection 
on the cyhnder is the affine center of one of its facets. Replace this facet of P with the cone 
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over Fn {x| ||x|| < 1} The part of the resulting polyhedral surface, that hes inside the sphere 
of unit radius, is indeed a PL-surface embedded in (in Klein's model). The surface if 
locally convex at every point and strictly convex at p. However, it is not the boundary of a 
convex body Notice that this surface is not simply connected. 

Consider a locally convex spiral, embedded in the y^-plane with two limiting sets: circle 
{p|a; — 0,y'^ + z'^ — 1} and the origin. That is this spiral coils around the origin and also 
around (from inside) the circle. Let M be the double cone over this spiral with apexes at 
(1,0,0) and (—1,0,0, intersected with the unit ball {x|||x|| < 1}. This non-convex surface 
is obviously simply connected, embedded, locally convex at every point, and strictly convex 
at all point of the spiral. 

A locally compact realization r of M is a realization such that for any compact subset C 
of X" C n r(M) is compact. The question remains: 

Problem 3 Is it true that any locally compact embedding of a simply connected surface in 
H" is convex? 

It remains an open question whether van Heijenoort's (1952) criterion works for embedded 
unbounded surfaces in H". We conjecture that it is, indeed, the case. The proof of the main 
theorem makes use of quite a number of technical propositions and lemmas. The proofs of 
these statements for M" by most part can be directly repeated for X", but in some situations 
extra care is needed. If the reader is referred to van Heijenoort's paper for the proof, it 
means that the original proof works without any changes. 

Notation: The caUigraphic font is used for sets in the abstract topological manifold. The 
regular mathematics font is used for the images of these sets in X", a space of constant 
curvature. The interior of a set S is denoted by {S), while the closure by [S]. The boundary 
of S is denoted by dS. Since this paper is best read together with van Heijenoort's (1952) 
paper, we would like to explain the differences between his and our notations, van Heijenoort 
denotes a subset in the abstract manifold M by S, while denoting its image in R" by S'; an 
interior of a set S in M" is denoted in his paper by S. 

The immersion r induces a metric on M by 

d{p,q) = GLBarcip,q)cM{\r{arc{p,q))\} 

where GLB stands for the greatest lower bound, and \r{arc{p,q))\ for the length of an arc 
joining p and q on M, which is the r-image of an arc joining these points on M. We will 
call this metric r-metric. 

Lemma 4 (van Heijenoort) Any two points of M can be connected by an arc of a finite 
length. Thus M is not only connected, but also arcwise connected. 

Lemma 5 (van Heijenoort) The metric topology defined by the r-metric is equivalent to the 
original topology on M. 

Lemma 6 (van Heijenoort) r{S) is closed in X" for any closed subset S o/M. 

Lemma 7 (van Heijenoort) If on a bounded (in r-metric) closed subset S G M, mapping r 
is one-to-one, then r is a homeomorphism between S and r{S). 
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The proofs of the last two lemmas have been omitted in van Heijenoort (1952), but they 
are well known in topology. 

Theorem 8 Let (n > 2) be a Euclidean, spherical, or hyperbolic space. Let M = (M, r) 
be an immersion of an {n — l)-manifold M in X", such that r(M) is hounded in X". Suppose 
that M = (M, r) satisfies the following conditions: 

1) M is complete with respect to the metric induced on M by the immersion r, 

2) 3yl is connected, 

2) M is locally convex at each point, 

4) M is strictly convex in at least one point, 

Then r is a homeomorphism from M onto the boundary of a compact convex body. 

Proof. Notice that our theorem for X" = H" directly follows from van Heijenoort's proof 
of the Euclidean case. Any immersion of M into H" can be regarded as an immersion into 
the interior of a unit ball with a hyperbolic metric, according to Klein's model. If conditions 
l)-4) are satisfied for the hyperbolic metric, they are satisfied for the Euclidean metric on 
this ball. Geodesies in Klein's model are straight line segments and, therefore, for a bounded 
closed surfaces in H", that satisfies the conditions of the theorem, the convexity follows from 
the Euclidean version of this theorem. 

The original Van Heijenoort's proof is based on the notion of convex part. A convex part 
of M, centered at a point of strict convexity o = r(o), o G M, is an open connected subset 
C of r(M) that contains o and such that: (1) dC = H (1 r(M), where H is a hyperplane in 
X**, not passing through o, (2) C lies on the boundary of a closed convex body Kc bounded 
by C and H. We call HflKc the fid of the convex part. Let Hq be a supporting hyperplane 
at o. We call the open half-space defined by Hq, where the convex part lies, the positive 
half-space and denote it by Hq. Wc call the r-preimage of a convex part C in M an abstract 
convex part, and denote it by C. In van Heijenoort's paper H is required to be parallel to 
the supporting hyperplane Hq of r(M) at o, but this is not essential. In fact, we just need 
a family of hypcrplanes such that: (1) they do not intersect in the positive half-space, (2) 
the intersections of these hyperplanes with the positive half-space form a partition of the 
positive half-space, (3) all these hypcrplanes are orthogonal to a line I, passing through o. 
Let us call such a family a fiber bundle {Hz}{i^Ho) of the positive half-space defined by I and 
Hq. (In the case of X" = M" it is a vector bundle.) In fact, it is not necessary to assume 
that I passes through o, but this assumption simplifies our proofs. Here z > denotes the 
distance, along the line I, between the hyperplane in this family and Hq. We will call z 
the height of H^. 

Proposition 9 A convex part exists. 

Proof, van Heijenoort's proof works for X", n > 2, without changes. ■ 

Denote by ( the least upper bound of the set of heights of the lids of convex parts centered 

at o and defined by some fixed fiber bundle {Hz}(^i^H^y Since r(JVt) is bounded, then ( < oo. 
Consider the union G of all convex parts, centered at o. We want to prove that this union 

is also a convex part. Let us depart for a short while (this paragraph) from the assumption 

that r(M) is bounded. G may only be unbounded in the hyperbolic and Euclidean cases. 

As shown by van Heijennort (1952), if X" = M" and C < oo, G must be bounded even when 
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r(M) is allowed to be unbounded. If X" = H" and ( < oo, G can be unbounded, and this is 
precisely the reason why van Heijenoort's theorem does not hold for unbounded surfaces in 
hyperbolic spaces. 

Since in this theorem r(M) is assumed to be bounded, G is bounded. Let us presume 
from now on that X" = §" (the case of X" = H" is considered in the beginning of the proof). 
dG belongs to the hyperplane and is equal to if^ fl M. dG bounds a closed bounded 
convex set D m H(^. Two mutually excluding cases are possible. 

Case 1: dimD < n — 1. Then, following the argument of van Heljenoort (Part 2: 
pages 239-230, Part 5: page 241, Part 3: II on page 231), we conclude that G U D is the 
homeomorphic pre- image of an (n— l)-sphere SU2) C M. Since M is connected, SUD = M, 
and M is a convex surface. 

Case 2: dimD = n — 1. The following lemma is a key part of the proof of the main 
theorem. Roughly speaking, it asserts that if the lid of a convex part is of co-dimension 1, 
then either this convex part is a subset of a bigger convex part, or this convex part, together 
with the lid, is homeomorphic to M via mapping r. 

Lemma 10 Suppose X" = Let C be a convex part centered at a point o and defined 
by a hyperplane from a fiber bundle {H;^}(^i ho)- Suppose B = dC is the boundary of an 
{n — 1)- dimensional closed convex set S in Hz- Either S is the r-image of an [n — l)-disk 
S m M and M = C U S, where C = r{C), or G is a proper subset of a larger convex part, 
defined by the same fiber bundle {Hz}[i^Ho)- 

Proof. Using a perturbation argument, we will prove this lemma by reducing the spher- 
ical case to the Euclidean one. Since S is {n — l)-dimensional and belongs to one of the hy- 
perplanes in the fiber bundle {Hz}[i^Ho): [conv G] fl Hq is either empty or (n — l)-dimensional. 
If it is non-empty, [conv C] fli^o must have a point other than o and its opposite. The closure 
of a convex set in X" is convex. Since [conv G] is convex, if it contains a point p of Hq other 
than o and its opposite, it contains some geodesic segment [op] lying in Hq. Since o is a 
point of strict convexity, there is a neighborhood of o on op all whose points, except for o, 
are not points of [convC], which contradicts to the choice of [op]. 

So, [convC] n Hq is definitely empty. Since, by Lemma El r(M) is arcwise connected, 
all of [convC], except for the point o, lies in the positive subspace. Therefore, there is a 
hyperplane H in §" such that [G] lies in an open half space if+ defined by H. We can regard 
Ei^ as a standard sphere in M'^"'"^. H defines a hyperplane in M'^"'"^. Consider an n-dimensional 
plane En in R"+^ parallel to this hyperplane and not passing through the origin. Central 
projection ri of M fl H^ on En obviously induces an immersion rir of a submanifold M' of 
M into En- 

This submanifold M' is defined as the maximal arcwise connected open subset of M 
such that (1) all points of this subset are mapped by r to and (2) it contains o. It is 
obviously a manifold. Let us prove that it exists. Consider the union of all open arcwise 
connected subsets that contain a. It is open and is acrwise connected, since it contains o. 
Let M' = (M',rir). 

The immersion rir obviously satisfies Conditions 2-4 of the main theorem |H| rir defines 
a metric on M'. Any Cauchy sequence on M' under this metric is also a Cauchy sequence 
on M under the metric induced by r. Therefore M' is complete and satisfies the conditions 
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of the main Theorem |H1 The central projection on En maps a spherical convex part of M on 
a Euclidean convex part of M'; it also maps the fiber bundle {Hz}(i^Ho) to a fiber bundle in 
the Euclidean n-plane En- van Heijenoort (1952) proved Lemma ITOl for the Euclidean case. 
Therefore, either M = C U S, or C is a proper subset of a larger convex part centered at o, 
and defined by the same fiber bundle {Hz}(i^Ho)- ■ 

The second alternative (C is a subset of a larger convex part) is obviously excluded, 
since C is the convex part corresponding to the height which is the least upper bound of all 
possible heights of convex parts. Therefore in Case 2 M is the boundary of a convex body 
which consists of a maximal convex part and a convex (n — l)-disk, lying in the hyperplane 



2 Locally convex PL-surfaces 

Theorem 11 Let r be a realization map from a compact connected manifold M of dimension 
n — 1 into X" (n > 2) such that M is complete with respect to the r-metric. Suppose that 
M = (M, r) is locally convex at all points. Then M = (M, r) is either strictly locally convex 
in at least one point, or is a spherical hyper-surface of the form S'^ ndC, where C is a convex 
cone in R"+^, whose face of the smallest dimension contains the origin (in particular, C may 
he a hyperplane in R^+^j. 

Proof. The proof of this rather long and technical theorem will be included in the full length 
paper of Rybnikov (200X). ■ 

Theorem 12 Let r he a realization map from a closed connected n-dimensional manifold M, 
with a regular CW- decomposition, in M" or (n > 2) such that on the closure of each cell C 
o/M map r is one-to-one and r{C) lies on a suhspace of dimension equal to dimC. Suppose 
that r(M) is strictly locally convex in at least one point. The surface r(M) is the houndary 
of a convex polyhedron if and only if each {n — 3) -face has a point with an M -neighhorhood 
which lies on the houndary of a convex n-dimensional set. 

Proof. M is locally convex at all points of its {n — 3)-cells. Suppose we have shown that 
M is locally convex at each A;-face, < A; < n — 3. Consider a. {k — l)-face F. Consider the 
intersection of Star(F) with a sufficiently small {n — k)-spheTe § centered at some point p of 
F and lying in a subspace complimentary to F. M is locally convex at F if and only if the 
hypersurface §nStar(F) on the sphere S is convex. Since M is locally convex at each k-face, 
§ n Star(F) this hypersurface is locally convex at each vertex. By Theorem ^2 S fl Star(F) 
is either the intersection of the boundary of a convex cone with § or has a point of strict 
convexity on the sphere S. In the latter case the spherical generalization of van Heijenoort's 
theorem implies that S fl Star(F) is convex. Thus M is convex at p and therefore at all 
points of F. 

This induction argument shows that M must be locally convex at all vertices. If is locally 
convex at all vertices, it is locally convex at all points. We assumed that M had a point of 
strict convexity. The metric induced by r is indeed complete. By van Heijenoort's theorem 
and Theorem M is the boundary of a convex polyhedron. ■ 
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3 New Algorithm for Checking Global 
Convexity of PL-surfaces 

Idea: check convexity for the star of each {n — 3)-cell of M. 

We present an algorithm checking convexity of PL-reahzations (in the sense outhned 
above) of a closed compact manifold M = (M, r). 

The main algorithm uses an auxiliary algorithm C-check. The input of this algorithm 
is a pair (T, T), where T is a one vertex tree with a cyclic orientation of edges and T is its 
rectilinear realization in 3-space. This pair can be thought of as a PL-realization of a plane 
fan (partition of the plane into cones with common origin) in 3-space. The output is 1, if 
this realization is convex, and otherwise. Obviously, this question is equivalent to verifying 
convexity of a plane polygon. For the plane of reference we choose a plane perpendicular to 
the sum of all unit vectors directed along the edges of the fan. The latter question can be 
resolved in time, linear in the number of edges of the tree (e.g. see Devillers et al (1998), 
Mehlhorn et al (1999)). 

Input and Preprocessing: The poset of faces of dimensions n — 3,n — 2, and n — 1 of M 
and the equations of the facets, OR the poset of faces of dimensions n — 3,n — 2,n — l, and 
of M and the positions of the vertices. We assume that we know the correspondence between 
the rank of a face in the poset and its dimension. There are mutual links between the facets 
(or vertices) of M in the poset and the records containing their realization information. All 
{n — 3)-faces of M are put into a stack Sn-3- There are mutual links between elements of 
this stack and corresponding elements of the face lattice of M. 

Output: YES, if r(M) is the boundary of a convex polyhedron, NO otherwise. 

1. while Sn-3 is not empty, pick an {n — 3)-face F from Sn-s] 

2. compute the projection of F, and of all [n — 2)-faces incident to F, 

onto an affine 3-plane complimentary to F; denote this projection by PStar(F); 

3. compute the cyclicly ordered one- vertex tree 7{F), whose edges 
are the (n — 2)-faces of PStar(F) and whose vertex is F; 

4. Apply to (PStar(F),7(F)) the algorithm C-check 

if C-check(PStar(F), T(-F)) = 1 then remove F from the stack Sn-3 

else Output :=N0; terminate 

endwhile 

5. Output:=YES 

Remark 13 The algorithm processes the stars of all {n — 3) -faces independently. On a 
parallelized computer the stars of all [n — 3) -faces can he processed in parallel. 

Proof of Correctness. The algorithm checks the local convexity of M at the stars of 
all (n — 3)-cells. M is compact and closed by Krein-Milman theorem (or Lemma ITT|) M has 
at least one strictly convex vertex. By Theorem local convexity at all vertices, together 
with the existence of at least one strictly convex vertex, is necessary and sufficient for M to 
be the boundary of a convex body. 
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Complexity estimates Denote by fk the number of /c-faces of M, and by fk,i - the 
number of incidences between A;-faces and /-faces in M. Step 1 is repeated at most fns 
times. Steps 2-4 take at most const /„_2,n-3 (Star (F)) arithmetic operations for each F, 
where const does not depend on F. Thus, steps 2-4, repeated for aU (n — 3)-faces of M, 
require 0(/„_2,n-3) operations. Therefore, the total number of operations for this algorithm 

is 0(/n-2,n-3)- 

Remark 14 The algorithm does not use all of the face lattice of M. 

Remcirk 15 The algorithm requires computing polynomial predicates only. The highest de- 
gree of algebraic predicates that the algorithm uses is d, which is optimal (see Devillers et al, 
1998). 

Prom a practical point of view, it makes sense to say that a surface M is almost convex, if 
it lies within a small Hausdorf distance from a convex surface S that bounds an n-dimensional 
convex set B. In this case, the measure of lines, that pass through interior points of B and 
intersect S in more than 2 points, will be small, as compared to the measure of all lines 
passing through interior points of B. These statements can be given a rigorous meaning in 
the language of integral geometry, also called "geometric probability" (see Klain and Rota 
1997). 

Remark 16 If there is a 3-dimensional coordinate subspace L o/M" such that all the sub- 
spaces spanned by {n — 3) -faces are complementary to L, the polyhedron can be projected on L 
and all computations can be done in 3-space. This reduces the degree of predicates from d to 
3. In such case the boolean complexity of the algorithm does not depend on the dimension at 
all. therefore, for sufficiently generic realizations the algorithm has degree 3 and complexity 
not depending on n. 

Remark 17 This algorithm can also be applied without changes to compact PL-surfaces in 
W orW. 
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